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We study the size dependence of the dynamical symmetry breaking in the four- 
dimensional Nambu-Jona-Lasinio model. We show that the presence of boundaries 
reduces the chiral breaking region, and this effect is strengthened for a larger num- 
I ber of compactified dimensions. A critical value for the length of the compactified 

■ dimensions exists, below which the dynamical symmetry breaking is not possible. 

Considering finite temperature and chemical potential, the chiral phase structure 
^ , for the system with compactified dimensions is obtained. A gradual decreasing of 

^ ■ the chiral breaking region with increasing of chemical potential is found. Also, at 



fixed chemical potential, the decreasing of the size of the system changes the order 
of the chiral phase transition. 
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1 Introduction 



The phase structure of the strongly interacting matter attracts a great deal of 
interest. Due to the intrincate mathematical structure of Quantum Cromody- 
namics (QCD), effective models that incorporate some of its properties have 
been largely employed. In this sense, four-fermion models, as the Nambu-Jona- 
Lasinio (NJL) model [1], are very useful for the investigation of dynamical 
symmetries when the system is under certain conditions, like finite temper- 
ature, finite chemical potential, gravitational field, etc. j2l[3lH] . In particular, 
an interesting aspect in the analysis of the phase transitions of four-fermion 
models is the effect of space compactification [5]|51l7f5f5|10fllfl2] . Boundary 
effects have also been considered for quark- meson models [T3fTi] . The general 
question is to estimate the relevance of the fluctuations due to finite-size effects 
in the phase diagram. With this purpose, different approaches have been used 
to study various aspects of these models, as the finite-size scaling analysis 
the multiple reflection expansion [TU], and the zeta-function method 



In this paper, we extend the techniques introduced in Ref. [TT] and inves- 
tigate finite-size effects on the dynamical symmetry breaking of the four- 
dimensional NJL model at finite temperature and chemical potential. This 
is done in the framework of zeta-function regularization and compactification 
methods [13]. This approach allows in a simple way to determine analytically 
the size-dependence of the effective potential and the gap equation. Then, 
phase diagrams at finite temperature and chemical potential, where the sym- 
metric and broken phases are separated by size-dependent critical lines, are 
obtained. Note that our approach leads to exact results within the mean-field 
(Hartree) theory without making use of any additional approximation. This 
allows to avoid, for example, difficulties appearing in calculating the density 
of states using the multiple reflection expansion, as pointed out in Ref. P^ . 



We organize the paper as follows. In Section II, we calculate the effective 
potential of the NJL model in the mean-field approximation, using the zeta- 
function method. The size-dependent gap equation is discussed in Section III, 
while the phase diagrams are shown and analyzed in Section IV. Finally, Sec- 
tion V presents some concluding remarks. 
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2 The formalism 



Our starting point is the massless version of the NJL model, described by the 
Lagrangian density, 

>C = gi^g + - E [{qa-^qf + [qil.a^qf] , (1) 

where q and q are the A^-component spinors, and the matrices a" are the 
generators of the group U{N), with — I/y/N. 

We perform the bosonization assuming that only one auxiliary field, that asso- 
ciated with the bilinear qa'^q, takes non- vanishing values. This auxiliary field, 
denoted here by a, plays the role of a dynamical fermion mass, such that when 
it has a non- vanishing value, the system is in the chiral broken phase. We shall 
work in the Euclidian space, performing the Wick rotation in the time coor- 
dinate. For generality, we consider the D-dimensional Euclidean space-time, 
restricting latter to the D — A case. Besides, taking a uniform, i.e. indepen- 
dent of coordinates, we obtain the effective potential up to one-loop order, at 
leading order in as 

where Ag/f is the effective action, V is the volume and 



Above, /iD is the dimension of the Dirac representation and A is a scale pa- 
rameter. 

To take into account finite-size effects on the phase structure of this model, 
the system is considered to have d (< D) compactificd dimensions. We denote 
the Euclidian coordinate vectors by a;^; and write xe = {y,z), where 



y 



(yi = Xe, yn = XE) , (^i = za = Xe) , (4) 



with the Zj component being defined in the interval [0, Lj] and n = D ~ d, 
corresponding to the topology M"'xS^"+^ x- • -xS^^ where S^^ is a circumference 
of radius Lj. Therefore, the compactification of the ^^-coordinates makes the 
/c^-components of the momentum A;^; to assume discrete values, 

K^j^{nj + Cj), (5) 
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where Uj = 0,1,2,... and Cj = | (j = 1,2,.. .,6?) for antiperiodic boundary 
conditions. Note that we may associate a given Lj with the inverse of tempera- 
ture (3 = 1/T, say Lj^ = /5, thus treating the system at finite temperature with 
d—1 compactified spatial dimensions. In this finite chemical potential, 

/i, can also be introduced through the rule = | — 

In the following, we use the zeta-function regularization method [15], rewriting 
Eq. ([3D as 

ho 



Ui(cr) 



where ({s) is given by 



2V 



C'(0) + lnA\(0) 



(6) 



ni,...^nfi=—oo 



K + K + ^' 



(7) 



where k1 = ijij + 0,)^, with V„ being the n-dimensional volume. The 

techniques of dimensional regularization can be used to perform the integra- 
tion over the non-compactified fc^-components, yielding 



C {s; {aj} , {cj}) 



(4vr)n/2 r (s) 



Yf s 



n 



;{a,},{c,} , (8) 



where F/^ (z/; {aj} , {cj}) is the generalized Epstein-zeta function, defined by 



17 (z.;{a,},{c,}) 



E cti ('^i + cif H h (rirf + Cdf + 



ni,...,n£j=— oo 



(9) 

with ttj = 47r^/L^. Note that F/^ is well-defined only for Re z/ > d/2, but 
it can be analytically continued, becoming a meromorphic function, into the 
whole complex i/-plane. 

As remarked in Refs. [T51ITT] . the analysis of the pole structure of the zeta- 
function implies that Eq. ([6]) must be written as 



Ui (<t; {aj} , {cj}) 



h 



D 



2Ki(47r)W2 



n 



T—r\ YT 



n 



^;{a,},{c,}), (10) 



for n odd, or 



Ui{a-{aj},{cj}) 



I'D 



-l)t 



2Vd(47r)"/2 



HI 
2 • 



Y? 



n 



;{aj},{cj} 



+ 17 



n 



;{aj},{cj} 



In A' 
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for n even, where 7 and V'('S) denote the Euler-Mascheroni constant (~ 0.5772) 
and the digamma function, respectively. Both expressions will be used in the 
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discussion of symmetry breaking in the four-dimensional NJL model with com- 
pactified dimensions. 



3 The gap equation 



The phase structure of the model is studied through the gap equation, obtained 
by minimizing the effective potential with respect to cr, 

= 0, (12) 

where m is the order parameter of the chiral phase transition, that is the 
dynamically generated fermion mass. We now discuss the size effects on the 
gap equation starting, for completeness and to set the free space parameters, 
by treating the model in absence of boundaries. 



d_ 



Ueff {a;{aj} ,{cj}) 



3.1 Non-compactified model 



Let us consider the NJL model in four dimensions without compactified di- 
mensions. This corresponds to d = 0, and so n = D = 4. This system has been 
treated in the literature (see Refs. [TTfTS] for reviews) using other methods. 
With the zeta-function approach, the renormalized effective potential is given 
by 

-Ueff{a,d = 0) = -—+ r 1-- A "T^TT-TnT^r 1 - - 



iV ' ' 2Gr (47r)^/2 V 2 J D(47r)^/2 V 2 

(13) 

where the renormalized coupling constant is defined by 

1 1 hoiD -1) / D\ n 9 

'-ril )X^-^. (14) 



Gr G (47r)^/2 V 2 

Note that Eq. ( |T3l) is valid for 2 < < 4; it is singular for D = 4, due to 
the pole of the gamma-function. In the limit D — ^ 4, we define e = 4 — D and 
obtain an analytically regularized effective potential, taking = 4, as 



1 Ueff (a) 1 6(t2 /I , . r 



-7 + ln47r + 



N A4 2GrX^ (47r)2A2 Ve 

+ (4;0^U-^ + ^^^^+2-^^p)^ 



this shows explicitly the singular behavior of Eq. (fT3l) as e — 0. 



Let us compare Eq. (fT5i) with the corresponding expression obtained by using 
the cut-off regularization [ T7p8] . 



1 U, 



cr 



N A4 



(47r)2A2 I P ~ 3 



+ 



a 



(47r)2A4 



^ A2 1 , 



(16) 

where A is the cut-off parameter which must be larger than A. Thus, the 
zeta-function and cut-off methods are equivalent through the correspondence 



- 7 In 47r + 1 



In 



A2 
A2 



(17) 



Hence, with the use of the correspondence in Eq. (ITTl) . the non-trivial solution 
of the gap equation derived from Eq. (fT5l) can be written as 



1 



1 

Go 



d 



UAa] 



Aw? , 

In — 



(47r)2A2 



mA^ c)cr 

U—'lli 

where we have defined the dimensionless coupling constant Gc = ^^Gr, and 
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(19) 



In Eq. ( ITS]) , it is possible to identify the constant Gq acting as a critical value; 
when Gc > Gq we have dynamically generated fermion mass. A value for Gq 
can be fixed by choosing values for the mass scale A and the cut-off A from 
phenomenological arguments; note that, for D = 4, the NJL model has to be 
considered as an effective model. 



3.2 Presence of boundaries 



To take into account temperature and finite-size effects, we must analyze the 
modified gap equation. 



(20) 



where m = m{{aj}, {cj}) is the boundary modified fermion mass. Then, ac- 
cordingly to the zeta-function approach, using Eqs. (fTOj) and (fTTj) . the modified 
gap equation, Eq. (!20|) . is written as 



a Go A2F,(47r)"/2 
while it has the form 



rfi-|)>1 



n 



- + 1; {aj},{cj}j , for d= 1,3; 

(21) 
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+ Yf (0; {oj} , {c,}) [in \^ - i - (1)] } , for <i = 2; 



and finally, 



1 
G 



1 



— = — - __ FP 



Yr (l;{a,},{c,}) 



for d = 4, 



(22) 



(23) 



where FP [1^4™ ] means the finite part of . 

Thus, taking the fermion mass approaching to zero in Eqs. (pTj) . (!22l) and (!23|) . 
we obtain the critical values of the coupling constant Gc with the corrections 
due to the presence of boundaries, for the cases d = 1,3, d = 2 and d = 4, 
respectively. In this context, YJ^ _^ ^ reduces to a homogeneous generalized 
Epstein zeta-function Y^. 



The next step is the construction of analytical continuations for Y^, which are 
written through a generalized recurrence formula [H] , 



r(.-i) 



Yd (i^; {%} , {cj}) = -^^^ — V ^^'^"^ V ~ 2' ^^^'^'^^ ' ^^■''^'^^ 

+ ^w^.(^-^;K},{9}), (24) 

where the symbol {aj^^} means that the parameter aa is excluded from the 
set {ttj}, and 

Wd {m Wj] , {cj}) = £ cos {27TndCd) {^=^ — \ Kr, ; 

(25) 



in the above equation X^_x = yJ2k~i '^k {nu + c^) and Ky{z) is the modified 
Bessel function of second kind (see Refs. [T91I20] for a discussion of the Cj ^ ai 
symmetry) . 



4 Phase structure 



4-.1 Finite-size effects on the critical coupling 



We now analyze the Lj-dependent critical curves for the phase diagram from 
the gap equation in the limit m — > 0. First, we study the compactification 



7 



of spatial coordinates at zero temperature. Thus, we take Eq. (12T|) in the 
cases d = 1 and d = 3, and Eq. (!22|) for the d = 2 case. After that, we 
use the recurrence formula, given by Eq. (^^, and perform the necessary 
manipulations. 



Considering the simplest situation in which all spatial coordinates are re- 
stricted to intervals of the same length and obey antiperiodic boundary con- 
ditions, i.e. Lj = L and cj = 1/2 for j = 1, . . . ,d, we obtain the following gap 
equations 

I _ 1 Ad 

a " " (Za)^' ^^^^ 

where 



Ai = -^ 0.16, 



As = 



oo oo 



A2 = A,-- j: 



"2 



TT 



7ii=— oo n2=l 



oo oo 



rii + 



n2 



K_i 27m2 



vr 



"i,n2=— oo n3=l 



1 

Til -\ 

2 



(27) 

0.22, 

(28) 



&^-2 



0.26. (29) 



To illustrate the results above, in Fig. 14. II we plot the critical coupling constant 
Gc, given by Eq. fl26l) . as a function of x = (LA)~^, with the coefficient A^ 
taking the values given by Eqs. fl27|) . (1251) and fl2^ : we also fix the value of Go — 
5.66, obtained from Eq. (pSD by choosing A ^ 1.25 GeV and A ^ 280 MeV 
The three cases considered correspond, respectively, to the system between two 
parallel planes a distance L apart; in the form of an infinite cylinder having a 
square transverse section of area L^; in the form of a cubic box of volume L^. 
For each situation the chiral breaking region corresponds to the region above 
the corresponding line. Since Gc increases as L decreases, we can infer that 
the finite-size effects require a stronger interaction to maintain the system in 
the chiral breaking region, as L is diminished. 

We also find from Eq. fl26|) that, for each d, a critical value Lc exists below 
which the chiral breaking region is completely suppressed; this critical value is 
obtained by setting the right hand side of Eq. fl2^ equal to zero. For Gq ~ 5.66, 
considering the values of A^ given by Eqs. fl27H2^ . we get = (LcA)~^ ~ 1.03 
for d = 1, Xc — 0.89 for d = 2 and Xc — 0.83 for d = 3. We see that the cubic 
box has the greatest value of which means that the largest number of 
compactified dimensions is the scenario which needs the strongest interaction 
to keep the system in the chiral breaking region. 

It is worth noting that the behavior of the critical coupling constant as a 
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1.0 H. 



0.8 - 




Fig. 1. The critical coupling constant Gc, Eq. (j26p . as a function of x = (LA)~^. The 
dashed, dotted and solid lines correspond to the cases of d = 1,2,3 compactified 
dimensions, respectively, with Aii given by Eqs. ([27ti29|) . We fixed Go ^ 5.66. 

function of L does not depend on the number of compactified dimensions, as 
it is expected from finite-size scaling arguments. 



4-2 System at finite temperature 



Now we consider the system with compactified spatial dimensions and at 
finite-temperature. Since the way of introducing finite-size effects is through 
a generalized Matsubara prescription, we can identify one of the compactified 
dimensions with the Euclidian (imaginary) time and take the compactification 
length as the inverse of the temperature, say Li = [3 = T~^. So, to analyze 
the (L, T)-dependent phase diagram, we must look at the critical equation 
obtained from Eqs. (fTgj) and (|20l). 



■ In — + =^ ^f/i {o- {a,} , {c,}) 



(47r)2A2 



0, 



(30) 



where ai = A-K^/Lf = Att'^/(3'^ and, for simplicity, we take 02 = • • • = = 
Att^/L^. Then, we look at Eq. ([22]) for = 2, Eq. ioi d = 3, and Eq. ([23]) for 
d = 4 considering, in all these cases, Li = (3. Next, we perform the necessary 
manipulations, in an analogous way as we have done in the previous subsection, 
and take the limit m — > 0. Thus, using Eqs. fl2^ and fl25]) . we obtain 



4m 



for d 



A2 4 ^ ^ 



\n2 



m 



ni- 



-00 712 = 1 
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(32) 



for d = 3; and finally 



4m 



E E( 

ni,"2,"3=— CO n4=l 



\n4 



/5 

xi^'i 27rn4— 




0, 



(33) 



in the case of c? = 4, where Ad, for d = 1,2, 3, are given by Eqs. (!27|) - (!29|) . 

In Fig. ??, we plot the phase diagrams corresponding to Eqs. (I3T1433I) in the 
{x, T)-plane, where x and T are the inverse of the compactification length and 
the temperature, respectively, measured in units of m, that is x = {Lm)~^ and 
T = (/3m) Each critical line separates the chiral breaking region, below the 
line, from the chiral restoration phase, above the line. 

In the limit x —>■ 0, corresponding to the system without spatial boundaries, 
the critical temperature is ~ 0.68 m. The size effects start to appear for 
X ~ 0.3. We see from the curves that the critical temperature decreases as the 
size of the system diminishes. Our results indicates that there is a minimal size 
of the system for the existence of a chiral broken phase. The critical values of 
the compactification lengths for suppression of the chiral breaking region are 
given by: Xc = {Lcm)~^ = 0.68 for d = 2, Xc = 0.58 for d = 3 and Xc = 0.54 
for d = 4. Note that, for the system between two parallel planes a distance L 
apart (the case d = 2), the critical line is symmetric by the change x ^ T, 
as expected. It can also be seen that the critical temperature has a faster 
decreasing with the size reduction for the system in the form of a cubic box. 
This fact is confirmed by the greatest value of the critical length (the length 
for which the critical temperature vanishes) in this case. 

In order to obtain a quantitative estimate for our results, we consider m as 
the constituent quark mass, m = 280 GeV according to Ref. [2T]. In this 
scenario, with the appropriate conversions, the critical temperature for the 
system without spatial boundaries is Tc ~ 0.68 m ^ 190 MeV and the size 
effects begin to be noted for L = {xm)~^ ^ 2.4 fm. For this choice of m, the 
critical values of the compactification length are Lc = (xcm)~^ ?s 1.1 fm for 
d = 2, Lr = 1.23 fm for d = 3 and = 1.32 fm for d = A. 
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4-3 System at finite temperature and at finite chemical potential 



The phase diagram can also be studied by taking into account the dependence 
on the chemical potential. As previously remarked, this is done by taking in the 
prescription for the momentum associated with the compactified Euclidian 
time, Li = (3 = 1/T and ci = | — Thus, considering the simplest situation, 
d = 1, which corresponds to the system in bulk (absence of boundaries), the 
use of Eqs. ( 1T8|) and (|2T|) leads to the following equation for the critical line, 

Am^ , 1 u2 

■ 111 ^ - 775? - TT^ = 0- (34) 



(47r)2 m2 6/52 27r2 

Now, if we consider the system with compactified spatial coordinates, the /i- 
dependence appears as additional factors in the last term of the left hand side 
of Eqs. (pTl) . ( |32l) and ( l33l) . which are obtained after the appropriate analytical 
continuation. 

So far we have considered only second-order phase transitions. However, the 
inclusion of a finite chemical potential can alter the nature of the chiral phase 
transition. To understand as this happens we proceed as follows. We expand 
the gap equation fl^Ul) for m near criticality, obtaining 

A{T, Lj, /i) + B{T, Lj,^)m^ + C{T, Lj, fx)Wi^ = 0, (35) 

where the coefficients A(T, Lj, fi), B{T,Lj,fi) and C{T,Lj,fi) are defined ac- 
cording to Eqs. (l2T]) - (!23l) . i.e. in terms of homogeneous generalized Epstein 
zeta-functions Yfi{i'), with different i/-exponents. The tricritical point is then 
determined by imposing the condition 

5(T,L„/i) = 0. (36) 



As an example, in Fig. 14.31 we plot the phase diagram in the {fi, T)-plane for 
the case without boundaries {d = 1). We identify the solid and dashed lines as 
representing, respectively, second- and first-order phase-transition lines. The 
dot locates the tricritical point, which occurs at Tc ~ 0.46 m ~ 128 MeV and 
fic ~ 0.90 m ^ 252 MeV. As expected, at high temperature or low chemi- 
cal potential, a second-order phase transition is suggested, while in the low- 
temperature or high-chemical potential limit the system experiments a first- 
order transition. 

Finally, in Fig. ?? we plot the phase diagram of the system with three com- 
pactified spatial dimensions {d = 4) in the {x, T)-plane, for some values of the 
chemical potential. The chiral breaking region is gradually diminished with in- 
creasing chemical potential. So, when n increases, lower values of temperature 
and smaller sizes are necessary to reach a phase-transition line. In addition. 
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Fig. 2. The phase diagram in the (^u, T)-plane for the case without boundaries 
{d = 1), with T and /i measured in units of m. Sohd and dashed hues represent 
the second-order and first-order phase transition hues, respectively, and the dot 
is the tricritical point. The chiral breaking region is below the line. Here we take 
A =4.46 EH. 

nonzero values of the chemical potential can alter the order of the phase tran- 
sition, as expected. Nevertheless, the decreasing of the size can also changes 
the order of the phase transition for the system at fixed /i. 



5 Concluding Remarks 



In this work, we have studied dynamical symmetries of a four-fermion model 
when the system is under certain conditions. We have investigated finite- 
size effects on the chiral phase structure of the four dimensional NJL model 
employing zeta-function regularization and compactification methods. In the 
mean field approximation, this allows to derive analytically expressions of the 
gap equation when (< 4) dimensions are compactified. 

First, the dependence of the critical coupling constant on the size of the sys- 
tem was analyzed, comparing the cases of one, two and three compactified 
spatial dimensions with the same compactification length L. The presence of 
boundaries implies that a stronger interaction is needed to maintain the sys- 
tem in the chiral broken phase; this effect is stronger for greater number of 
compactified dimensions. It is shown that there is a minimal size of the sys- 
tem below which there is no chiral breaking region; the corresponding critical 
values of the compactification lengths, for ci = 1, 2, 3, are determined. 
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The chiral phase structure at finite temperature was investigated by taking 
one of the compactification lengths, that associated with the compactificd 
Euclidian time, equal to the inverse of temperature. The results suggest that, 
as the size of the system diminished, finite-size effects start to appear for 
a given value of the compactification length, with the temperature rapidly 
decreasing as 1/L is further increased; this behavior is more accentuated for 
the case with the greater number of compactificd spatial dimensions, giving 
the largest decreasing of the chiral breaking region. 

The study is concluded with the analysis of the dependence of the chiral phase 
transition on the chemical potential. For the bulk system, i.e. the system with- 
out spatial boundaries, we get a simple expression for the phase-transition line, 
showing that the critical temperature decreases as the chemical potential in- 
creases. The nature of the phase transition changes from second- to first-order 
at specific values of /i and T which locates the tricritical point. For the system 
with compactificd spatial dimensions, we find that the chiral breaking region 
decreases as the chemical potential is increased; that is, when /i increases, 
lower values of temperature and smaller sizes are necessary to keep the sys- 
tem in chiral broken phase. We show that, for a fixed value of ji, the decrease 
of the size leads to a change in the order of the transition. 

Possible extensions of this work are the investigation of properties of fermion- 
fermion condensates under finite-size conditions, as well as, the study of the 
system under the infiuence of an external magnetic field. 
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